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Abstract 

We propose a method to evaluate hadronic matrix elements of QCD gluon oper- 
ators in the instanton vacuum. We construct the ground state of the interacting 
instanton ensemble for non-zero ??-angle using a variational principle. A method to 
study the i9-dependence of observables on the lattice is suggested. We then derive 
the effective fermion action, which allows to calculate hadronic correlation functions 
in a 1/Ac-expansion (Nambu-Jona-Lasinio type effective fermion theory). Gluon 
operators are systematically represented as effective fermion operators. Physical 
matrix elements are obtained after integrating the correlation functions over fluctu- 
ations of the numbers of instantons. The influence of the fermion determinant on 
the topological susceptibility is taken into account. Our effective description gives 
matrix elements fully consistent with the trace and U{1)a_ anomalies. The approach 
allows to consistently evaluate the nucleon matrix elements of various gluon and 
mixed quark-gluon operators in a chiral soliton picture of the nucleon. 
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1 Introduction 



It is now widely recognized that instantons play an important, if not crucial, role in the 
dynamics of strong interactions. After their discovery twenty years ago their relevance 
to hadron phenomenology was explored in a number of pioneering papers 0, 0, while 
a quantitative theory of the instanton medium was still lacking. A treatment of the 
interacting instanton ensemble by means of the Feynman variational principle showed how 
the instanton medium stabilizes itself at rather low densities 0. In particular, the ratio 
of the average size of the instantons to the average separation between nearest neighbors, 
p/R, was computed as ~ 1/3, the value suggested from phenomenological considerations 
1^. The instanton vacuum was also studied in direct lattice experiments, by the so-called 
cooling procedure^] J^, ||, ^ |10|, |Tl|, [1^. It was observed that instantons and antiinstantons 
(J's and /'s for short) are the only non-perturbative gluon configurations surviving after a 
sufficient smearing of the quantum gluon fluctuations. The measured properties of the // 
medium, such as the ratio p/ R, appear to be close to those computed from the variational 
principle 1,0,0. 

The success of the instanton vacuum lies in its explanation of all phenomena related to 
the dynamical breaking of chiral symmetry in QCD. The mechanism of chiral symmetry 
breaking is the delocalization of the would-be fermion zero modes of individual instantons 



T3| , p^ , leading to a finite fermion spectral density at eigenvalue zero, proportional to the 
chiral condensate, (ipip) [jl^. Equivalently, one observes that quarks propagating in the 
instanton medium acquire a dynamical momentum-dependent mass []T1|, |TB|. A massless 
pion appears as a collective excitation. 

It should be noted that instantons do not lead to confinement . However, it has been 
realized in many ways that it is chiral symmetry breaking, and not confinement, which 
determines the basic characteristics of light hadrons (except , probably, for highly excited 
states) 1^, O, [1^. This point of view has recently received direct support from lattice 



measurements of correlation functions in the cooled vacuum [|TT|. Cooling the quantum 
fiuctuations eliminates not only the perturbative one-gluon exchange, but also the linear 
confining potential. Nevertheless, the correlation functions of various hadronic currents 
in the cooled vacuum, where only J's and J's are left, appear to be quite similar to those 
of the true "hot" vacuum. 

Detailed numerical studies of correlation functions of meson and baryon currents in the 
instanton medium by Shuryak, Verbaarschot and Schafer have shown impressive agree- 
ment with phenomenology [jl9|, Correlation functions in the instanton vacuum can 
also be computed analytically [|14|, |17|. By integrating first over the instanton coordi- 



nates, one derives an effective fermion action of the form of a Nambu-Jona-Lasinio model 
T7[]. It exhibits a many-fermion interaction with a specific spin-fiavor structure, first 



suggested by 't Hooft p[. The effective fermion theory can be studied systematically in a 
l/Ai'c-expansion. In this formulation, chiral symmetry is broken by the instanton-induced 
many-fermionic interactions. Correlation functions of meson currents can be computed 
as fermion loop diagrams with a momentum cutoff defined by the average instanton size, 
p"^ ~ 600 MeV. The study of correlation functions of baryon currents at large leads 
to a picture of baryons as chiral solitons — N^. valence quarks moving in a self-consistent 



-•^For a more complete list of references, see [||. 
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meson field |]2T|. This approach to low-energy baryon structure has been worked out 
extensively and gives a good account of the static properties of non-strange and strange 



baryons [^]. Furthermore, by calculating correlation functions of baryon currents with 
various weak and electromagnetic current operators, baryon form factors have been com- 
puted [^, which are found to be in excellent agreement with experiment. 

In the description of the deep-inelastic structure of the nucleon, one faces the task 
of evaluating nucleon matrix elements of certain operators containing gluon fields |2^]. 



There are essentially two types of such gluonic operators. Operators of twist 2 are related 
to the moments of the nucleon structure functions. Operators of higher twist, such as 
F^^, F^uFf^i^, or ipFfj,i,j^ip arise in the description of power corrections. Their evaluation 
is a challenge for any model of the QCD vacuum capable of describing the nucleon. Of 
particular interest is the nucleon matrix element of FF, the subject of the so-called proton 
spin problem^ [|^. The nucleon matrix elements have to be extracted from correlation 
functions of nucleon currents with the respective gluon operators. A related problem, 
which also requires evaluation of matrix elements of operators involving fermion and gluon 
fields, is the estimation of vacuum condensates of "mixed" operators like V'-F^vCr^^?/', which 
give corrections in QCD sum rules [p^ . 

The object of this paper is the study of hadronic matrix elements of gluon operators 
in the instanton vacuum. Our framework is the grand canonical ensemble of interacting 
/'s and /'s, to which fermions are coupled mainly through the zero modes. The fluctua- 
tions of the number of /'s (A^+) and /'s (A^-), are essential to realize the renormalization 
properties of QCD, such as the trace anomaly and the f/(l)^-anomaly. Our strategy 
will be as follows. We first consider the instanton ensemble in gluodynamics and derive 
the distribution of the sizes and numbers of pseudoparticles from a variational principle 
1^. With this instanton size distribution, we determine the effective fermion action in 
quenched approximation, for a fixed number of instantons, N±. Correlation functions 
are then evaluated in a two-step procedure. In the first step, we average over instanton 
coordinates and the fermion field for fixed N±. For correlation functions of operators 
consisting of fermion fields (such as meson or baryon currents), this average can be per- 
formed directly with the help of the effective fermion action. Moreover, as we shall show, 
even correlation functions with gluon operators can be reduced to averages in the effective 
fermion theory: in the course of integrating over instanton coordinates, the QCD gluon 
operator is replaced by an effective many-fermion operator. In the final step, we then 
pass from canonical to grand canonical averages, by integrating over fluctuations of the 
number of J's and J's, N±. The probability for these fluctuations is given by the instanton 
partition function including the fermion determinant. By this procedure, we can calculate 
correlation functions of nucleon currents with gluon operators, and extract the nucleon 
matrix elements by taking appropriate large-time limits. 

Our intention here is twofold. First, we want to set up a framework for concrete calcu- 
lations. All correlation functions are systematically reduced to averages in a Nambu-Jona- 
Lasinio type effective fermion theory, which can be evaluated using established techniques 
|2^ . Second, we want to demonstrate that the fundamental renormalization properties 



of QCD are preserved in this approach, even at the level of hadronic matrix elements. 



2 



For a recent lattice estimate of the nucleon matrix element of the topological charge, see EE 
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For this reason, we shall focus on the operators and F^yF^y. In QCD, the nucleon 
matrix element of is uniquely fixed by the trace anomaly, while the matrix element 
of F^yF^y is related to the isoscalar axial charge by the U{1)a anomaly. In the instan- 
ton formulation, the values of these operators are related, respectively, to + and 
A^_i_ — A^_. We shall show that the nucleon matrix elements calculated according to our 
prescription are obtained in agreement with the QCD anomalies. The proper inclusion of 
non-topological and topological fluctuations of the numbers of instantons will be crucial 
for this property. 

Matrix elements of gluon operators depend in general on the renormalization point. 
When computed in the instanton vacuum, they are obtained at the scale set by the 
instanton background, 1/p. In general, the results must be subjected to QCD evolution, 
to make contact with experimental data at higher momentum scales. This problem will be 
addressed elsewhere. Here, we concentrate on gluon operators related to QCD anomalies, 
about which exact statements can be made in QCD at the hadronic scale. 

The plan of this paper is as follows. In section 2, we consider the grand canonical 
ensemble of instantons in gluodynamics. The distribution of the total number of instan- 
tons, + A^_, is derived from the QCD trace anomaly. We discuss the dependence of 
physical quantities on the -i^-angle and suggest a way to measure this dependence on the 
lattice. We then perform a variational estimate of the grand canonical partition function 
at finite ^9-angle, allowing for a general form of average pseudoparticle interaction with 
different strengths between same and opposite-kind instantons. In section 3, we include 
fermions in the partition function. We derive the effective fermion action in zero mode 
approximation, for an ensemble with a fixed number of instantons, N±, but allowing for 

7^ A^_. We outline the bosonization of the 't Hooft interaction in the case of more than 
one fermion flavor {Nf > 1). We then discuss the suppression of the topological suscepti- 
bility due to light fermions in QCD on the basis of the U{1)a and chiral Ward identities. 
We show that the fermion determinant obtained in the instanton vacuum for ^ 
precisely describes this suppression. In section 4, we consider correlation functions of 
hadronic currents with gluon operators. We show how fixed- correlation functions can 
be represented as averages over the effective fermion theory, in which the gluon operator 
is replaced by an effective fermion operator. We formulate the rules how to "translate" 
an arbitrary gluon operator into an effective fermion operator, and demonstrate their 
consistency. The correlation functions of the grand canonical ensemble are then obtained 
by averaging over fluctuations of + and — A^_. It is shown that the basic 
renormalization properties of QCD are preserved in this framework: we verify the low- 
energy theorem of scale invariance for correlation functions with the operator / d^x F^^, 

and the realization of the U{1)a anomaly for the operator / d^x F^yF^y. In section 5, we 
then derive the nucleon matrix elements of and F^yF^y in our approach. The nucleon 
matrix element of F^j, is obtained consistent with the trace anomaly. The matrix elements 

of F^yF^y reduces to the isosinglet axial coupling constant, g^^ i calculated within the 
effective fermion theory. This allows us to consistently interpret the result of a former 
chiral soliton calculation of q'a^ as an estimate of the nucleon matrix element of F^yF^y 
in the large-iVc limit ||2^. Conclusions and an outlook are presented in section 6. 
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2 Statistical mechanics of instantons 



2.1 The grand partition function of instantons 

The essence of the instanton vacuum is that the number of "pseudoparticles" is not fixed, 
hence one is deahng with a grand canonical ensemble. Let us denote the number of J's 
and /'s, respectively, by A^+, A^_. The grand partition function is generically given as 

Z{fx,'d) = J2 exp{{fx + i^)N+)expi{fx-i^)N_)ZN^, (2.1) 

N+,N- 

l N++N- 

Zn± = ^ ,^ , n J d'^zjdpidUi do{pi) exp{-Uint)- (2.2) 

Here, the integrals go over the collective coordinates of instantons: their center coordi- 
nates, zj, sizes. Pi, and orientations, given by SU{Nc) unitary matrices in the adjoint 
representation, Uj; dUj means the Haar measure normalized to unity. It is useful to in- 
troduce complex chemical potentials for /'s and /'s, p ± i-d, where i) is the usual vacuum 
angle. The instanton interaction potential, f/int, depends on the separation between pseu- 
doparticles, zj — zj, their sizes, pi,pj, and their relative orientation, UjUj. By do{p) we 
denote the one-instanton weight. In the one-loop approximation it is given by 0, |2^ 

doip) = %/3(Me,02^=exp[-/?(p)], (2.3) 



p^ 

where P{p) is the one-loop inverse charge. 



= ^ = Mog(J^). ^, = H^,^. (2,4) 



Note that f3 in the pre-exponential factor of eq.( p.3|) starts to "run" only at the two-loop 
level, hence its argument is taken at the ultra-violet cut-off, M^ut- The coefficient C^^ 
depends on the renormalization scheme; in the Pauli-Villars scheme it is 

4 60 e-^-^^^'^ 

= vr2(iV.-l)!(iV.-2)!- ^^"^^ 
The relation between the QCD scale parameters. A, in different schemes is 

Apv = 1.09 A^js = 31.32 Aiattice. (2.6) 
If the scheme is changed, one has to change the coefficient according to Cat C'^ = 

In the two-loop approximation, the instanton weight is given by |2^ 



doip) = %/3(p)^^^exp 



(2.7) 



where Pu{p) is the inverse charge to two-loop accuracy. 



Mp) = m + |log^' b' = Iat^. (2.8) 
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The partition function, eq.( |2.1| ), is actually normalized to the perturbative partition func- 
tion (without instantons). Therefore, when the four-dimensional volume, V, goes to 
infinity, which we always assume, the logarithm of Z gives the ground state (vacuum) 
energy density, 6*44, with that of perturbation theory subtracted. 



exp [-V (6*44 - (6'44)pert)] • 



(2.9) 



Here, ^^i, is the energy-momentum tensor. If the vacuum state is isotropic, one has 
6*44 = \^9^^. Using the trace anomaly. 
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Pig 



^9 



4 flU — 



32n 



2 M^^' 



(2.10) 



where F^^ = F^u^^u is the gluon field strength squared and (3{g^) is the Gell-Mann-Low 
function (not to be confused with the inverse charge /3(p)), 



dg\M, 



cut ) 



d log Mcut 



-h- 



Stt^ 



b' 9' 
2 (87r2)^ 



(2.11) 



one gets 



exp V 



b 1 



^2 \ 



4 32vr2 r ^-/np^ 



(2.12) 



Here, (-^^j,)np is the gluon field vacuum expectation value, which is due to non-perturbative 
fluctuations, i.e. the gluon condensate It is important that this is a renormalization 
group-invariant quantity^, meaning that its dependence on the ultraviolet cutoff, M^ut, 
and the bare charge given at this cutoff, (7^ (Mcut), is such that it is actually cutoff- 
independent. 



327r2 



72 \ 

tj.u/ 



np 



Mcut exp 



3^ (A/cut) dg^ 



1 4 



Pig' 



c'Mc^ut exp 



327r2 



(2.13) 



(Here, c, c' are constants independent of Mcut-) By definition of the gluodynamics partition 
function. 



^glue 



the l.h.s. of eq. ( |2.13| ) is equal to 



72 \ 

327r2 V^'''/ 



np 



V(M, 



cut J 



d ^ I ) 



(2.14) 



1 1 



d log Zgiue 



VM^d[l/g^{M,^,)\ 



(2.15) 



^To be more precise, the RNG-invariant quantity is (0^^), see eq.( ^.lO ). However, if the coupHng 
constant at the ultra-violet cutoff scale, f;^(Afcut)j is small enough, it is sufficient to use the beta function 
to one-loop order. 
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Applying the same differentiation to eq.( |2.13| ) once more, one gets a low-energy theorem 
11, 



d'^ log Zgiu 



1 



(87r2)2d[l/^2(M,„0 



(327r2)2 
4 1 



d^xF, 



6 327r2 



I np 



(327r 



2^2 



/ np 



(2.16) 



np 



If the bare coupling, g^iMcux)^ is not chosen small enough, there are obvious corrections 
to this formula, following from the higher-order terms in the beta function, eq. (|2.11|) . 

This low-energy theorem has important consequence for instantons: it predicts the 
dispersion of the number of pseudoparticles in the grand canonical ensemble 0. Assuming 
the instanton ensemble to be sufficiently dilute, and taking into account that the one- 
instanton action is 



1— inst 



one can rewrite the low-energy theorem, eq.( |2.16| ), as 



12 



■{N), 



(2.17) 



(2.18) 



UN, 

where = A^+ -|- A^_ is the total number of pseudoparticles. Thus, it follows directly 
from the renormalization properties of the Yang-Mills theory that the dispersion of the 
number of instantons is less than for a free gas of instantons, for which one would get a 
Poisson distribution, (A^^) — (A^)^ = (A^)- At Ac — > oo, the dispersion becomes zero. 

One concludes that the interaction of instantons with each other is crucial to support 
the necessary renormalization properties of the underlying theory: any cutoff of the in- 
tegrals over instanton sizes "by itself", such as a cutoff due to a possible infrared fixed 
point of the beta function, as suggested recently by Shuryak leads to the Poisson 
distribution and hence to the violation of the low energy theorem, eq.( |2.18| ). 

By further differentiating logZgiue with respect to the bare coupling, g'^{Mcut), one 
can easily generalize eq.( p^.l8D to higher moments of the distribution. For example. 



{N^) - 3{N^){N) + 2(A^)^ 



4X2 



(2.19) 



etc.. Summarizing these formulas, one concludes that the distribution of the number of 
pseudoparticles should, for large (A^), be given by 



P{N) oc exp 




- 1 



(2.20) 



A convenient way to generate all moments of the number distribution is to use the chemical 
potential, /i, introduced in eq.( |2.1| ). Renormalizability predicts that the grand partition 
function, eq.(|2.1|), as a function of /x, should behave as 



Z(/x) 



exp 



(2.21) 



Differentiating log Z{fj,) with respect to fi and setting /i 
in accordance with the low-energy theorems. 



0, one obtains all the moments 
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2.2 7?— angle dependence 

The statements made in the previous section are independent of the -(^-angle and should 
thus hold at any This means that in eqs.( p.20| , |2.21|) only the quantity (A^) may depend 
on ■}). Hence the partition function can be written in the form 



Z{fi, ■&) = exp 



/(O) = 1, (2.22) 



where (A^)o is the average total number of pseudoparticles at -i? = 0. The function f{i)) 
can not be determined from general considerations. Nevertheless, comparing eq. (|2.22| ) 
and eq.( p.l| ), one realizes that eq.( |2.22| ) imposes a non-trivial restriction on the instanton 
interactions. The factorization of the /x- and "i^-dependence would be a property of non- 
interacting particles, but in this case the factors 6/4 would be replaced by unity. 

Information on /(i?) can be extracted by differentiating logZ with respect to The 
second derivative is proportional to the topological susceptibility of the vacuum. 



= {{N^~N^r) = -l{N)o^y. (2.23) 

The dependence of the QCD observables on the angle is of fundamental importance. 
However, until now no direct lattice measurements with non-zero i} have been performed. 
The reason is that the inclusion of a non-zero 'd-teim in lattice simulations is ambiguous 
and, moreover, makes the integration measure complex. The above discussion suggests 
an unambigious method to study the -(^-dependence of observables on the lattice. Indeed, 
cooling down a gluon configuration, one finds — and that rather early — the topological 
charge of the configuration, Qt = A = — N_, which is the canonical conjugate to i). 
If -P(A) is the fraction of configurations with given A, the probability of a configuration 
with given -(9 is given by the Fourier transform, 

Z{^) = 5]P(A) exp{i^A). (2.24) 

A 

Measuring a QCD observable, O, such as the string tension etc., for a set of configura- 
tions, and knowing the topological charge of those configurations, one can establish the 
dependence of the vacuum average of the observable on -i?, 

0{^) = J20{A)P{A) exp{i^A). (2.25) 

A 

The topological charge of a configuration, A = — A^_, can rather easily be estab- 
lished after cooling. The total number of /'s and /'s, = -|- A^_, is a less well-defined 
quantity, as it tends to decrease during the lattice cooling. The question arises: When to 
stop the cooling procedure? When does one hit the true instanton vacuum in the lattice 
simulation? The answer is provided by the dispersion of A^, eq.( p.l8 ), or, more generally. 



by eqs.(|2.20"|, |2.21|), which differ from the Poisson distribution by factors of 6/4 = 11 Nc/12. 
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When that distribution in the total number of pseudoparticles is achieved, one can measure 
various observables hke the condensates, correlation functions, etc.. The most non-trivial 
check, however, that one is dealing with the true instanton vacuum, is the abovementioned 
factorization of the /x- and '(9-dependence in the logarithm of the partition function. Let 
the probability of a configuration with A^_|_ J's and /'s be P(A^+,iV_). If the origi- 
nal gluon configurations are cooled sufficiently, but not over-cooled, the grand partition 
function should have the form of eq.( |2.22| ). 



N± 



exp 



^-{N)o fm exp(^ 



(2.26) 



Once such non-trivial behavior in /i and is established, one can read off the gluon 
condensate as 

^(f-i) = (2.27) 

3271^ \ /np Y 

The above statements are about all one can say on the dispersion of pseudoparticles 
on general grounds. To obtain more information on the instanton vacuum, one has to go 
into dynamics, i.e., take into account the instanton-instanton interaction in a quantitative 
way. In the next section, we estimate the ground state of the grand canonical partition 
function, eq.(|2.1|), using a variational principle. 



2.3 Variational estimate of the grand partition function 

In this section, we show that actually very limited information on the instanton inter- 
actions is needed to make a variational estimate of the basic properties of the instanton 
ensemble. We shall neglect the many-body forces, an approximation which is justified a 
posteriori, if the medium turns out to be sufficiently dilute. 

An intuitively clear description is achieved by introducing the actual size distributions 
of instantons in the medium, which do not necessarily coincide with the weight of an 
isolated instanton, do{p), of eqs.(pl3|, p7|) . Moreover, the results of section 2.1 show that 
they necessarily do not coincide with (io(p) — otherwise the renormalization properties 
of the theory would be violated. Following ref. 0, we now derive the "best" distribution 
functions, d±{p), from a variational principle. This approach preserves all the general 
requirements discussed in sections 2.1 and 2.2. 

As an intermediate step, it will be convenient to work with the canonical ensemble. 
We first derive the best d±{p) for a given number of J's and J's, and leave the summation 
over N± for the end. Let us write the partition function for given N±, eq. (|2.2|) , in the 
form 

-, N++N. N++N_ 

Zn^ = J^^J n di:eM-mi m = - E log(io(p/) + t/int(0- 

(2.28) 
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Here, C,i = {zi,pj, Uj) denote the set of collective coordinates of the Jth pseudoparticle, 
and ^ those of all + A^-_ pseudoparticles. 

For the purpose of a variational estimate, we introduce a "non-interacting" partition 
function, 

. iV+ N. . N++N. 

= iV+!FTyn^e/c^+(p/)n^e/c^4Pj) = ^^^^ n d^expl-E^iOl 

EliO = -Elog^+(P/)-Elog^-(P/)- (2.29) 

Here, d±{p) are the effective size distributions in the ensemble of interacting /'s and 
J's. Note that these distributions, naturally, depend on N±. The partition function, 
eq.( |2.29| ), is immediately expressed through fugacities, C±, 

Zi = nJn1(^(-^)''''(^^-)''-^ C± = / dpdM- (2.30) 

To find the best d±{p) for a given interaction potential, Ui^t, we use the Feynman vari- 
ational principle, which is based on the convexity of the exponential function. It states 
that 



= Ziexp[-{E - El)] > Ziexp[-{E - El)]. (2.31) 
Here the averaging is performed with the modified partition function, Zi, 



A{0 = ^^^nivT/n^e/^+(P/)n^e/^-(P/)^(0. (2.32) 



ZiN+lNJ.J'-j^ ^ 
In our case we have 



E-Ei = E/ J dplog 



Ce -J do{p) 



1 ^ N N r 
+ E 7^7^/«2 4,(pi)4.(p2)t/int(6,6). (2.33) 



The best effective one-particle distributions, d±{p), are those which maximize the r.h.s. 
of the inequality eq.( |2.31| ), i.e., those which reproduce best the true partition function. 
We see from eq. (|2.33|) , that in this approach the optimum distribution depends only on 
averages of the interaction potential over separations and orientations. There are two 
basic quantities, the average interaction between instantons of the same and of opposite 
kind; different average interaction between II and // are prohibited by CP invariance. 
Let us denote the basic averages by 

J dUi dU2 d^zi d^Z2 ullS''\zi - pi, P2, UlUl) = Vf3p\pl^a, (2.34) 
j dUi dU2 d^zi d^Z2 Uli^"\zi - Z2, Pi, P2, UlUl) = V(3pipl^s- (2.35) 
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Some comments are in order here. The factors V and (3 on the r.h.s. of eqs.( p.34| , |2.35| ) 
are trivial consequences of the fact that the potential depends only on the separation, 
and that it is proportional to the inverse charge, f3. Less trivial is, which argument one 
should one choose for the running f3. Following the reasoning of 0, we are inclined to 
write P{{R)), where {R) is the average separation, to be determined below. Fortunately, 
this question is not very important quantitatively, as the dependence on the argument 
is only logarithmic (see also the discussion of the numerical values below). The product 
of pIp2 arises because of dimensions, and because the interaction has to vanish at zero 
sizes. Finally, we have introduced here two phenomenological constants for the average 
interaction of pseudoparticles of the same and opposite kind, 7^ and 7^. It should be 
mentioned that for the sum ansatz one has |^ 



7a 



27 . 
7s = -r^' 



Nr 



m - 1 



(2.36) 



However, here we would like to take a more liberal stand and allow arbitrary 7s, 7a. 
Nevertheless, the behavior of 7s,7a at large as l/Nc, shown by eq. (|2.36| ), is probably 
a general statement — it arises due to averaging over relative orientations in eqs.( |2.3^ , 
PI). 

Substituting the definitions eqs.( |2.34 |2.35| ) into eqs.( |2.30| , |2.31| , |2.33| ), we get 



> 



^exp(^J2N, log(V^Ce)-log 



doip) 



p 



X exp I -— ^^^l^^2PlPile,e 



(2.37) 



Here, we have introduced the average sizes, 

1 



Pi 



dp d±{p)p^ 



(2.38) 



Varying eq. (|2.37|) with respect to d±{p) we find the best effective size distribution, which 
takes into account the interactions with the medium, 



d±{p) 



Q:_i 



C±rfo(p) exp(-p2a±). 



V 



a. 



V 



(2.39) 

(7aP^iV+ + %P^Ar_). (2.40) 



The coefficients C± cancel in eq.( p.37D and can be both put to unity. 

We see that the one-instanton size distributions are multiplied by gaussian cutoff 
functions, whose slope depends on the density of particles. 
eqs.( p.30| , ^.3^ ) one obtains the fugacities, (±, 



Substituting eq.( |2.39| ) into 
and the average sizes of the instantons, p^. 
In this way, we get the best variational approximation to the fixed- A^-t partition function. 



> 



1 



/5 



exp|^iV,log(\/C) + ^ E N,.N,2PlPlle,e 



(2.41) 
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Taking the one-loop expression for do{p), eq.(|2.3|), one finds 

. o . 1 



pI 



j dpdo{p) 

aiogCi 



2 

6-4 



/3 = /5(M,,t). 



(2.42) 



(2.43) 



da± a± ' 2 
We get a system of equations relating the average sizes of J's and J's to their densities, 



(2.44) 



At = N_ = N/2, which will be the case at "i? = 0, when the vacuum is CP symmetric, 
one has the packing fraction. 





( 




V 


Pi 


[isPl^ 




^' 




( —N+ 




V 






+ Ispl^) 





{A) 



2 N 

V 



2 N 
V 



13 Is + la 



(2.45) 



To get the variational estimate of the fixed- partition function, Zjv±, one has to 
substitute p± from eq. (|2.44|) into eq. (p.41|) . Our aim is, however, the grand partition 
function, 2'(/i, ■(9), which is obtained from Z^^ according to eq. (|2.1|) . The summation over 
N± in eq. (|2.1| ) can, of course, be done by the steepest descent method. The saddle point 
values of A^4- and are in general complex and conjugate to each other, and both are 
proportional to the four-dimensional volume, V — * oo. At "i? = 0, the saddle point values 
of N± coincide, and are real. After some straightforward algebra we obtain the logarithm 
of the grand partition function, 

logZ(/i,t9) = ^(iV)o exp(^)/(^9). (2.46) 

Here, (A^)o, the average total number of J's and J's at -(9 = and = 0, is given by 



P 



(2.47) 



Note that we have obtained the grand partition function exactly in the form predicted in 
section 2 by considering the renormalization properties of the Yang-Mills theory. Here, 
we get a concrete formula for the density of instantons in terms of the Aqcd parameter, 
and a concrete form of f{'&)- The function f{i)) in eq. (|2.46| ) is given in an indirect way. 
One has to first solve the equation determining the phase of the saddle point of N±, 
xlj = \og{N+/N.)/2i, 



ip + — - — arctan ■ 



^9, 



cos-?/' + rA/l — sin ijj 



Ja 
Is'' 



(2.48) 
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substitute its solution to determine the saddle point value of the partition function, and 
then read off f{i^), 



r cos ip + yl — r"^ sin i/j 
1 + r 



4-i) 
b 



(2.49) 



This function depends on the ratio, r, of the average interactions of same and opposite 
kind pseudoparticles, eqs.( |2.34| , p.35| ). 

It should be noted that not all values of 7s,a lead to a stable ensemble. First, 7^ + 
7a must be positive — otherwise the system would collapse. This can be seen from 
eqs.( p.44| , p.47|) . Second, 7^ itself must be positive. If the interaction between same- 
kind pseudoparticles were on the average attractive, the system would prefer to break 
up into two subsystems, with only J's in one part of the universe and only J's in the 
other, with an apparent maximum violation of CP symmetry in both parts. Third, the 
topological susceptibility of the vacuum, eq.( |2.23| ), must be positive. Determining the 
second derivative of f{i)) from eqs. (|2.48| , |2.49| ), we find 



{Qt) 



h-r{h- 4) 



(2.50) 



Hence we must require r < b/{b — 4), meaning that the effective repulsion of opposite-kind 
pseudoparticles can not be much stronger than that of same-kind. 

In several simple cases, the function f{-t}) can be found analytically. For example, at 



1, i.e., 7a = 7s, which is the case of the sum ansatz of 0, we have 



(2.51) 



and the topological susceptibihty is 

(Qt) = {N)o- (2.52) 
In the case r = (no effective repulsion between opposite-kind instantons), we have 



f{^)r=0 = COS—. 

At r —1, which is the edge of the stability region, we get 



(2.53) 



2t9 



2(6-2) 



COS ■ 



6-2 



(2.54) 



By measuring the T^-dependence of the partition function by lattice cooling, as described 
in section 2, one should be able to "experimentally" determine the ratio r. This quantity 
is one of the basic characteristics of the instanton vacuum. 

To end this section, we present numerical values of the basic characteristics of the 
instanton vacuum, which were calculated using the two-loop instanton density, eq.( p.7| ). 
One first determines the average instanton sizes from eqs.( p.38| , p.44j) . The average in- 



stanton density then follows from maximizing eq.( |2.41| ). Since the argument of the inverse 
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Qj 


/ N y/^ 

WAV 


1 / y1 

V v) 


1- 1 / /I -1 


I 


2.5 


.68 


.37 


5.8 




3 


.66 


.35 


7.5 




3.5 


.64 


.34 


9.0 




4 


.63 


.33 


10.2 




5 


.62 


.31 


12.4 




7 


.60 


.30 


15.8 


II 


2.5 


.74 


.42 


6.8 




3 


.72 


.40 


8.5 




3.5 


.71 


.39 


10.0 




4 


.70 


.38 


11.3 




5 


.68 


.36 


13.5 




7 


.66 


.34 


16.9 



Table 1: The dimensionless instanton density, {N /V A.'^y/^ , the packing fraction, 
(l? N/Vf/^, and the effective inverse charge, [3, for various values of a parameter a, 
which defines the argument of the inverse charge according to (3 = h\og a{N /V K^Y^^ . 

The two sets correspond to different values of the constant \{'^a + 7s) determining the 
average instanton interaction, cf. eqs.( p.34| , |2.35|) . Set I is obtained with the sum ansatz, 
eq. (|2.36|) , set II with |(7a + 7s) equal to half the value of the sum ansatz. 



charge, /?, is not uniquely determined from first principles, we set (3 = 6 log a^N/VA^Y^"^ 
with a a dimensionless parameter, and explore different choices of a. Results for the instan 
ton density (i.e., the gluon condensate), {N/VA^Y^*^ the packing fraction, (p^ N/VY^'^, 
and the effective inverse charge are shown in table 1. Given there are the results for two 
different values of the constants defining the average instanton interactions, 'ya,s- Here, 
= iV_, so that the results depend only on the sum, 7a + 7s- From table 1 one sees that 
the instanton density and packing fraction depend only very weakly on the argument of 
the inverse coupling, a, although the values of (3 vary considerably. The values obtained 
for the density and packing fraction are rather close to those of 0. 



3 Including fermions 

3.1 The effective fermion action 

So far we have discussed the partition function of instantons in pure gluodynamics. When 
fermions are included {Nf flavors), the coefficient of the beta function, eq. (|2.4| ), is changed 
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to 



b = —N^--Nf. 
3 3^ 



(3.i: 



The fixed-A^± partition function with fermions, normahzed to the perturbative partition 
function, is given as 



Y fermions 



Det(m, Meut) 



1 



, N++N. 



NJ.N_ 



-J n c?6c?o(P/)exp(-f/i„t)Det(m,Me,t), (3.2) 



det(iy (0 + im) det{i^ + iM^^t) 



det {i^ + im) det (iy (0 + iM^ 



(3.3) 



Here, m = diag(mi, . . . , m^r^) is the bare quark mass matrix, M^nt the Pauh-Villars 
regulator mass. Furthermore, y (^) denotes the fermion Dirac operator in the background 
of the iV±-instanton configuration. 



y(0 = <^~tAix;0, 



(3.4) 



(Here, A± are the fields of individual J's and I's, in singular gauge 0.) The eigenvalues 
of the Dirac operator contributing to the determinant can can be divided into "low" and 
"high" frequencies by introducing a splitting mass parameter. Mi, 



Det (m,M 



cut , 



Det(m,Mi)Det(Mi,M, 



cut , 



(3.5) 



The value of Mi will be chosen of order 1/p. The essence of the instanton description of 
chiral symmetry breaking is to treat the two factors in eq. (|3.5|) in different ways 0, [1^ . 
The high-frequency part can approximately be factorized in one-instanton contributions, 
and is therefore included in the one-instanton weight of eq. ( |3.2|) . The low-frequency 
part, which is dominated by the fermionic zero modes associated with the instantons, 
must be averaged in the background of all N± instantons simultaneously. This leads to 
the delocalization of the zero modes, which is the mechanism of chiral symmetry breaking. 

When the fermion determinant is included in the variational description of the inter- 
acting instanton partition function, it modifies the effective instanton size distribution in 
the medium (in addition to the change in the /3-function, eq.( ^.l| )). This modification is 
of the order Nf/N^ To lowest order in this parameter (quenched approximation), 

one gets back the unmodified size distribution of gluodynamics, eq.( ^.39| ). In the follow- 
ing we use the gluodynamics size distribution in one-loop approximation to average the 
low-momentum fermion determinant. With the sum ansatz, eq. (p.36|) , the effective size 
distributions for /'s and J's coincide, and eq.( |2.39| ), with eq. ( |2.3| ), becomes 



d±{p) 



const X ^ exp 



6-4 p2 



Pi 



(3.6) 
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Since b = 0{Nc), the normalized size distribution reduces to a (5-function in the hmit 



Nc oo, 



J dp±d±{p)^ d{p) S{p-p) {b oo), p = \ffp. (3.7) 

For simphcity, we shall use the 5-function distribution in the following calculations^, i.e., 
we replace all instanton sizes by p. 

The variational approximation to the fixed- A^-t partition function with fermions can 
then be written as 



^fermions ^ Dct^^ , Detjv± = / H -^^f// Det (m. Mi ) . (3.8) 

Here, Zj^^ is the variational partition function of gluodynamics. We remark that correc- 
tions to this factorized form of the partition function can in principle be calculated using 
the variational principle. 

In the ground state described by the partition function, eq.(|3.8|), the chiral symmetry 



of the fermions is spontaneously broken O, 14 1. A useful concept is the effective fermion 



action [17|. It is obtained by performing in eq.( |3.8|) the integral over the instanton coordi- 
nates, and is valid in the vicinity of the saddle point of eq.( |3.8D . In terms of the effective 
action, the fixed-iV-|- partition function represented as 

D^ATi = C j Vil^^V^ exp (-S,^[ij\i;]) . (3.9) 

(Here, C is a normalization constant.) The effective action can be used to calculate 
averages of operators consisting of fermion fields (see section 4). 

To explicitly construct the effective fermion action, we need to know the fermion 
propagator in the background of the A^-j. -instanton configuration, eq. (|3.4|) . For a single 
instanton, the fermion propagator is singular in the chiral limit due to the zero mode. It 
may be approximated as the sum of the free propagator and the explicit contribution of 
the zero mode [|I4|, W%, 



+ ^my' ^ - ^^("'^^(^-))^^^^'^^(^-)). (3.10) 

J l-inst ^ ^ ^ im ^ ^ 



Here, $±(x; ^/(/)) is the wave function of the fermion zero mode in the background of one 
/(/). This interpolating formula should be accurate both at small momenta {p ^ 1/p), 
where the zero mode is dominant, and at large momenta [p ^ l/p), where the propagator 
reduces to the free one. Eq. (|3.1Cl| ) is equivalent to using an approximate fermion action 



m 



exp(-5^(^")[^t^^]) 0, explj2jd'xtp}i^^f\Yi{imf-Vl'^'\'ip},tPf]),{3.n) 



f / /=i 



^This approximation is only a practical simplification, and does not imply that we take b to infinity 
elsewhere. Corrections for finite-width size distributions can readily be included. 
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Vi^'\^},^Pf] = |ci^x(^}(x)z^<l>±(a;;6(/)))/A(*i(2/;6(J))^^^/(y))- (3-12) 

The fermion Green function calculated with this action coincides with eq. (|3.10|) . Given 
the diluteness of the instanton medium, one may thus approximate the fermion action in 
the background of an A^-|--instanton configuration by the product 



exp (^—S[ilj\ilj]j oc exp / d'^x tp^j^ic^ipf 



(3.13) 



/ 

N+ Nf 



X 



I /=1 



n n {^^f - viii^h v^/]) 1 1 n n (^^/ - vii^h ^/] 



^N- Nf 



I /=1 



With the fermion propagator approximated by eq. (|3.13|) , the average over the config- 
urations of the independent instantons, eg. (|3.8| ), reduces to a product of averages over 
individual /'s and /'s and can be carried out explicitly. Since the splitting mass, Mi, 
is of the order of the effective momentum cutoff in the zero mode wave function, 1/p, 
the normalizing determinant of eq. (|3.3|) in the zero mode approximation may simply be 
replaced by Mi to the power of the number of zero modes. 



'++N-) 



(3.14) 



Its only role in the following will be to make the partition function dimensionless. The 
averaged fermion determinant, eq.(p.8|), becomes 



N++N. 



^Nf{N++N^) 



J V^^Vijj exp J d^x^/j^i^ijjf] W++ W^- . (3.15) 



Here, W± denote the one-instanton averages 



V 



Nf 



f 



(3.16) 



Performing the average over color to leading order in l/N^, eq. (|3.16|) can be expressed asQ 



— detJ±(2;) + ^irrif 



Nf 

E 

/=i 



2-2\Nf-l 



,Nf 





fd^z^ 




[ nJ J 





iJ±iz) + 



mNr 



J -^det;j±(z) 

(3.17) 

+ 0(m2). (3.18) 



^Explicit expressions for the zero mode vertex, eq.( [3.12 ) and formulas for integrals over the color 
orientation matrices can be found in O] . 



18 



Here the determinants are over flavor indices, and detj denotes the minor in which the 
/-th row and column are omitted. The currents, J±{z), are color singlets and Nf x Nf- 
matrices in flavor, 



J+(x, 



fg 



/ exp(-^(A: - l).x) F{k)F{l) ^t(fc)i(i ± ^,)^^(/), (3.19) 



where F{k) is a form factor related to the wave function of the zero mode in momentum 
space. 



Fik) 



1 t^o 
t 00 



(3.20) 



Ik p. 



The first term of eq. (|3.17| ) has the form of the 't Hooft interaction . The form eq. (|3.18| ) 
is equal to eq. (|3.17|) to leading order in m. 

To determine the effective fermion action, we have to find the saddle point of the 
fermion integral, eq. (|3.15| ). We introduce an integral representation for the powers of the 
one-instanton averages. 



Det 



N± 



exp log 



'47rV 



\+V \N^Mi 



/Nf 

J Vip^Vip exp J d^x ij2'^h^^f + ^+ 



iJ+{x) + 



iV+ + (+ 
niNr 



(3.21) 



+ + 



We now determine the common saddle point of eq. ( |3.21| ) in the variable \± and the 
fermion field, in the thermodynamic limit, N± 00, V — >■ 00, N±/V fixed. 

The case Nf = 1. Let us first consider the case of one quark flavor, Nf = 1. In this 
case, the exponent in the fermion integral of eq. (|3.21|) is bilinear in the fermion field. 



''%\k)\{l±r.)mF\k) + ^,(3.22) 



(27r 



and the fermion integral can be carried out exactly. 



Det7v± 
W^i(A+,A_) 



"2^ 



d\^ 



27r 



exp(-iyi(A+,A_)) 



iV+log 



N^ 



2^2 ~ 



- iV, + (+ 



mVNc , , 



— Tr log 



^ -1^ [A+(l+75) + A-(l-75)]F^(A:) 
Ijt — im 



(3.23) 



(Here, Tr includes the trace over momentum space, and color and spinor indices.) Since 
the fluctuations of A = N^ — N_ are small relative to = A^+ + A^_ in the thermodynamic 
limit, it is sufficient to evaluate the effective action to leading order in A. We set 



A 



± 



M(l ±5), 



(3.24) 
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at the saddle point, with 6 = 0(A). The values of M and S are determined by minimizing 
the effective potential, eq. (|3.23|) . This leads to the condition 



{2TxYM^F\k) + k'^ 



N 



mMVNc 
27r2p2 



(3.25) 



This equation is identical to the self-consistency condition of |]14[. Note that M tends to 
a finite value in the thermodynamic limit. Parametrically, 



Furthermore, 



M oc 



S 



N 



1/2 



P- 



mMVNr 



at the saddle point. The effective fermion action, eq. (|3.9|) , is thus determined as 



Nf=l 



(2vr) 



-ip^ik) Ijt -iMF\k){l + 6-f5) ipik). 



(3.26) 



(3.27) 



(3.28) 



It describes fermions with a dynamical momentum-dependent mass, which is proportional 
to the square of the wave function of the instanton zero mode. For A 7^ 0, i.e., for an 
unequal number of /'s and J's, the fermion action shows a parity-violating mass term. 
This term diverges in the chiral limit; we have retained only the infrared-singular part 
and dropped all finite terms. This infrared singularity is canceled when averaging over A 
in the grand canonical ensemble and does not affect infrared-stable physical quantities 
(see section 4.2). 

The case Nj > 1. In the case of more than one quark flavor, Nf > 1, the instanton- 
fermion vertex, eq. (p.l8|) , describes a many-fermionic interaction {2Nf fields). In order 
to perform the integral over the fermion field, we have to introduce auxiliary boson fields 
to linearize the exponent of eq.( |3.21| ). This can be done in the large- Ai'c limit, using the 
formula 



exp(A det [iJ]) 



dM. exp 



-{Nf-l)\ ''f-'{detMYi'' +iii{MJ) 



(3.29) 



which holds in saddle point approximation. Here, is a hermitean Nj x A'^^-matrix 
variable. With the help of eq. ( |3.29| ), eq.( |3.21| ) is represented as 



Det 



N± 



dX^ 



dX. 



exp A^_|_ log 



Nf 



X+V \NcMi 



-) (3.30) 



VM+VM-exp / d^x 



iNf~l)X+'''"\detM+ix))''f-' - (+ 

Nr. 



+ 



47]-2p2 



X 



tr [m{M+{x) + M-{x))]j 
Vip^V^ exp J d^x \ + [M+{x)J+{x) + M-{x)J_{x)] 
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One can now integrate over the fermion field. Furthermore, one may integrate over the 
auxihary variables, A±, in saddle point approximation, and represent the partition function 



as 



DetjVi = JvM+VM-exp{-W[M+,M-]) 



(3.31) 



For constant meson fields, A4±{x) = M.±, the saddle point condition for A± is 

Nf-l 



det A^±. 



(3.32) 



Assuming the vacuum meson field to be diagonal in flavor, A4± = diag(A^±i, . . . J^±Nf] 
the effective potential is obtained as 



W{M+,M-] 



iV+log 



M+fV [n,M, 



-iV, + (+ 



(3.33) 



f 



Tr log 



^ - \i [M+f{l + 75) + ■M-/(l - 75)] F\k) 
Ijt — im 



To leading order in Nc, the fermion determinant is simply the product of the determinants 
for the different fermion flavors, and the effective potential is the sum of the potentials 
for one flavor, eq. (|3.23|) . Setting 



M 



(3.34) 



the minimization of eq. (|3.34|) thus leads to the same conditions as for Nj = 1, eqs. (|3.25| , 
|3.27| ), for each flavor. Substituting this vacuum value of the meson field in eq.( |3.32| ), and 
returning to the representation eq.( p.21|) , we obtain the effective fermion action in the 
chiral limit, 



Ses[ip\^] = -Ud'xJ2^}t^i^f + {1 + 5)Y+ + {1-5)y\, 



5 



Nf-l 



iMff J d^x det J±(x), 



/ 



27r2p2A 

MVNr 



(3.35) 

(3.36) 
(3.37) 



This general formula includes as a special case the result for Nf = 1, eq.( p.28|) . The 
dynamical fermion mass is the same as for Nf = 1. (Again, we have dropped 0{m) terms 
in the CP-symmetric part of the action, so that Mi = . . 



Mjv, = M.) 
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3.2 Fermions and the topological susceptibility 

In deriving the effective fermion action for given N±, we fiave used tfie instanton size 
distribution of pure gluodynamics, eq. (|3.6| ), an approximation justified by tlie formal 
parameter Nf/Nc. Indeed, the modifications of the size distribution due to the low- 
energy fermions are only quantitative. However, when considering the grand partition 
function, eq. (|2.1| ), there is one important aspect in which fermions qualitatively alter the 
picture. The presence of dynamical light fermions modifies the probability of topological 
fluctuations, A = — N_ 7^ 0, in a principal way: (A^) must vanish in the chiral 
limit. This remarkable property is a consequence of the ?7(l)^-anomaly of QCD. In the 
following, we derive the topological susceptibility from the QCD anomaly equation and 
the chiral Ward identities. We then show, that the instanton vacuum leads to an identical 
result in the chiral limit. 

The anomalous divergence of the isosinglet axial current operator is given by 

Nf _ N ~ - 

E^'(Mm75^/) = -/^FF + 2tY.mf^f^,^f. (3.38) 
/ / 

For the time-ordered product of two topological charge densities, one obtains the Ward 
identity 



2«Em/yc/^xT(0|FF(x), ^n5^/(0)|0) = ^ j d^xT{Q\FF{x),FF{S))\Q) . 

(3.39) 



/ 



To determine the topological susceptibility, we saturate the left-hand side of eq. ( |3.39| ) 
with pseudoscalar meson states. In general, there are N'j pseudo-scalar bosons related 
to the breaking of U{Nf) x U{Nf), of which Nj — 1 become massless in the chiral limit, 
(Goldstone bosons), while one singlet state remains massive. The couplings of these states 
to the non-singlet and singlet axial currents can be written as 

9^(01 ^t"7M75V'(0)|vrA(g)) = F^MIX% (a = 0, . . . , iVj - 1). (3.40) 

Here, the flavor generators, t", are normalized according to 



tr 



(3.41) 
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A unitary mixing matrix, X^, takes into account that, for m 7^ 0, mass and flavor 
eigenstates are in general not identical. (The index A runs over the N'j "nonet" meson 
states.) We now need the matrix elements of the operators on the left-hand side of 
eq. (|3.39|) between the vacuum and one-meson states. The topological charge density. 



FF, is a flavor singlet, and its matrix element can be parametrized as 

iFF(O)lvr^) ^ -^Xl (3.42) 
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where k is a constant, to be determined below. The matrix element of the pseudoscalar 
density, 'ipt"''^^ilj, is obtained from the non-singlet Ward identity, 

-2 /rf4xr(0|V^{r,m}+75V^(x), V^tSV^(O)lO) = (0|7/;{nt^}+V^(0)|0). (3.43) 



Inserting the pseudoscalar mass states, tta, and using the definition of the mixing matrix, 
eg. (13.401) , one finds 



(0|Vit'^75V'|vr^) 



t — v\ A . 



(3.44) 



Here, 



(3.45) 



Finally, the mixing matrix, X^, is determined by taking the matrix element of the anomaly 
equation, eq. (|3.38|) , and the non-anomalous divergence of the non-singlet axial current. 



d^ii'Wt''^) = #{t^m}+75V^ (a = l,...,iVj-l), 
which leads to an equation 



-tr 



(3.46) 



(3.47) 



Inserting now singlet and non-singlet states into eq. (p.39|) , using eqs. (|3.42| , |3.44| ), one 
obtains 



167r2 



rf^a;T(0|FF(x),FF(0)|0) 



2N 



f 



A Ml 



(3.48) 



Solving eq. ( ^.47| ) for X^, one obtains after a straightforward calculation 



V 



F2 



'Nf 



(3.49) 



In the chiral limit, m — 0, this becomes 

Ml = _ 

V 



'Nf 



-1 



(3.50) 



This replaces the gluodynamics formula, eq. (|2.23| ), in the presence of light fermions. The 
topological susceptibility is proportional to the harmonic average of the quark masses, 
which means that it vanishes if (at least) one of the fermion flavors becomes massless. 
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It is interesting to consider eg. ( |3.49| ) also in the "quenched" hmit, i.e., to switch off 
the fermions before letting the quark masses go to zero. (Formally, this corresponds to 
the limit Nj 0, with S"^/^ = 0{Nf) and (ipip) = 0(1).) In this case, eq.( |3.49| ) reduces 
to 



^ /('5t)quenched 



(3.51) 



Taking now the limit m — 0, the mixing between singlet and non-singlet states disap- 
pears, and one obtains from eq.(|3.47|) 



singlet ■ 



K = m:„„,„,. (3.52) 

Eqs.( p.51 , |3.52|) are the well-known Witten-Veneziano formula for the topological suscep- 
tibility in quenched QCD 



In the instanton language, the vanishing of the topological susceptibility for m — 
can be understood as being due to the "unbalanced" fermionic zero modes in the fermion 
determinant at A 7^ 0. To see this explicitly, we calculate the fermion determinant in zero 
mode approximation, eq. (|3.8| ), as a function of A. Evaluating the integral representation, 
eqs.( p.31| , |3.33| ), at the saddle point defined by eqs.( |3.25| , |3.34| , |3.37|) , we obtain, to leading 
order in 1/m, 



Det 



oc exp 



N,MV 



(3.53) 



On the other hand, the quark condensate, calculated from eq. (|3.33|) at A = 0, is 



Nf 



Minkowski 



Euclidean 



N. 



E 



1 d 



V dm 



f 



log Det 



N± 



A=0 



27r2p2- 
(3.54) 



One thus sees that the A-distribution described by the fermion determinant of eq.( 
precisely realizes the theorem eq.(p.50|). 



NcMV 
27r2p2 



-1 



' Nf 



m 



f 



(3.55) 



Note that the A-distribution corresponding to the full variational partition function, 
eq. (|3^ ), including the gluonic part. 



P(A) oc Z 



fermions 
N± 



Z^^DetN^ oc exp 



A2 



2N . 



exp 



2V{iJij) y ^ 



(3.56) 



is dominated by the one from the fermion determinant, since 



N 
V 



'Nf 



» -(^^) E 



m. 



(3.57) 
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in the chiral limit. 

The variational partition function of independent instantons, eg . (13 .81 ), including the 
fermion determinant, thus correctly describes the dispersion of topological vacuum fluc- 
tuations, which in QCD is governed by the U{1)a anomaly. Contrary to the case of the 
trace anomaly, which is realized already in the quenched approximation (see section 2.1), 
the inclusion of the fermion determinant is essential here. In section 4, we shall see that 
the independent instanton partition function with fermions, eq.( |3.8| ), can be used also to 
evaluate matrix elements of the topological charge in hadronic states consistently with 
the U{1)a anomaly. 



4 Evaluating gluon operators 
4.1 Effective operators for gluons 

To extract information about hadrons, one has to evaluate correlation functions of hadro- 
nic currents, together with some operators consisting of fermion and gauge fields. The 
prescription for calculating vacuum averages of functions of the gauge field and the fermion 
field is defined by the grand canonical instanton partition function with fermions, with 
the fixed- partition function given by the variational approximation, eq. (p.8|) . The 
grand canonical averaging is performed in two steps. First, for a fixed number of J's and 
/'s, one integrates over the collective coordinates of the instantons and over the fermion 
field, taking into account the dynamically generated quark mass. In the final step, one 
then averages over the number of J's and J's according to their dispersion in the medium. 
In this section, we investigate fixed-A^± (canonical) averages of operators, specifically of 
such involving gluon fields. The averaging over the number of J's and J's will be discussed 
in the following section. 

Let us consider first averages of operators involving only fermion fields. For example, 
the baryon correlation function is defined as the vacuum average of two baryon currents, 
consisting of Nc fermion fields coupled to a color singlet, 

M^) = ^^'^^•■■""^r.,;,...,,^^,^C/.(^)---0..(^)' (4-1) 

(r is a spin-flavor matrix) |21|. For a fixed number of instantons, N±, the normalized 
average over the instanton coordinates and the fermion fields can be represented as an 
integral with the effective fermion action, eq. ( |3.35D , 

'Jn{Xi)JI!{x2)] 



fixcd-Af± 



= Detj,lJv4j^V^ Mxi)Jl{x2) J n -YdUiexp(-S[ilj\^]) (4.2) 
ViI)'^ViIj JN{xi)jlf{x2) exp (-S'ef^[?/'^ 



j Vij^Vij exp(-Seff[V^^V^] 



= {JN{Xl)JN{x2))eS- (4.3) 



Here and in the following, we use the symbol (. . .)efr to denote the normalized average of 
a function of the fermion field within the effective fermion theory, as defined by eq.([4.3|). 
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The fermion integral can be evaluated to the order in l/N^ to which one knows the effective 
action. 

A new situation arises when calculating averages of operators containing gluon fields. 
Consider the nucleon matrix element of a gluon operator, e.g., the local gauge- 

invariant operators F^yF^j_i, or F^yF^^., or functions thereof. The nucleon matrix element 
is obtained from the correlation function of T[A\ with nucleon currents, eq.( |4.1| ), 

(4.4) 



fixod-Ar± 

N++N- 



Here, the function T[A\ is evaluated with the gauge field of the A^-|--instanton configu- 
ration, eq.( p.4|) , and averaged over the instanton coordinates together with the fermion 
action in the instanton background, eq.( p.l3|) . We want to represent eq.( |4.4| ) as an average 



within the effective fermion theory of the form of eq.(4.3): 



Jn{xi)A{^2)HA])^^^^_^^ (4.5) 

pv^tp^jiv(xi)4(x2) exp(-^eff[v^^^]) 



P^tp^exp(-S'eff[^^^]) 
J^(a:l)J^(x2)"^"[V^^V^]' 



' eff 

Here, ^^J^''^[%p\ip\ is a (generally non-local) fermion operator, which represents the gluon 
operator, in the effective fermion theory. We shall call "J-"" [t/'^, ^/^j the effective 

operator to T[A\ and denote it by the symbol of the gluon operator put in quotation marks. 
As with the effective fermion action, eq. ( |3.35| ), the definition eq.( |4.5| ) is understood in the 
sense of the 1/iVc-expansion, i.e., the effective operator is defined only in the vicinity of 
the saddle point of the effective fermion theory. 

There are several advantages in performing the averaging in the order of eq.( [4.5| ), 
integrating over instanton coordinates first. One obtains a simple and transparent repre- 
sentation of correlation functions with gluon operators. Furthermore, once the effective 
operator for a given gluon operator has been determined, it can be inserted in correla- 
tion functions of many different hadronic currents. The resulting fermionic correlation 



functions can be evaluated in the chiral soliton picture of the nucleon |21 



To construct the effective operator, we have to perform in eq.(|4.4|) the integral over 



the coordinates of the N± instantons, and cast the remaining fermion integral in the form 
of an integral with the effective fermion action, eq.( p.35|) . Since the effective operator is 



defined only in combination with the effective fermion action, it is important that eq. ([4.4|) 
is evaluated consistently with the integral defining the effective fermion action, eq.( |3.9|) , 
using the same zero-mode approximation for the fermion propagator, eq. ( |3.13| ). 



For an iV±-instanton configuration, eq. (|3.4|) , the function of the gauge field, ^[^], can 



be decomposed in one-instanton and multi-instanton contributions, 

N+ Af_ 

= E-^[^+(e/)] + E-^[^-(e/)] + 2-inst. + 3-inst. + ... (4.6) 
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When averaged over configurations of independent instantons, tlie multi-instanton terms 
are proportional to additional powers of the packing fraction of the instanton medium, 
p'^N/V ^ relative to the one-instanton terms. As the instanton medium is dilute, these 
contributions are numerically small. We retain only the one-instanton contributions in 
eq.( ^.4| ). It will be seen that this is consistent with the approximations made in deriving 
the effective fermion theory. 

In the following, let us consider for simplicity eg. ( [4. 41 ) without the baryon currents, 
i.e., the vacuum average of ^[^]. With the zero mode approximation for the fermion 
propagator, eq. (|3.13|) , we obtain in analogy to eq.( |3.15|) 



fixcd-Af± 

.N++N. ^i^^ /N++N. 



Det-^|p^tp^ J ll'^dUil^ E J'[A±mjexp{-S[^p\^P]) (4.7) 
D^^^ J V^P^ViP exp 1^ j d'^xijy^ijf 



X 



Again we have used the fact that the instantons are independent. Here, W± are the free 
one-instanton averages, eq. (p.l6|) , while WV± denote the one-instanton averages with an 



insertion of the function evaluated with the gauge field of /(/) 

"^'^dUj^j) ^[A±(6(/))J 

/ 



74 _ J^/ _ 

W^^[^P\^|J] ^ (-l)^^7^?^^^^(^")-^[^±^^^(^"))] n^i^'^tV^}'^/] + 0(m).(4.9) 



Here, J^[A] is averaged over the collective coordinates of one /(/) together with the zero 
mode fermion vertex, eq.( |3.12D , which leads to a coupling of the gluon operator to the 



fermion fields. In the important case that J^[A] is a color-singlet function of the gauge 
field, it is independent of the color orientation of the instanton. The average over orien- 



tations is thus the same as in the free instanton-fermion vertex, eq.( 3.17 ), and leads to a 
fermion vertex similar to the 't Hooft interaction, 

Yr±[ij\ij] = ( — ] (zM)^/yrfV)-^[^±K/)'f^/(/) = l'P)] detJ±(^).(4.11) 



(The prefactors have been chosen in analogy to the definition of the vertices, Y±, eq. (|3.36| ). 



If J^[A] is not a color singlet, which may be the case if it is part of a mixed fermion-gluon 
operator the average over color orientations in eq.( 4.11| ) results in fermion vertices with 
non-singlet color structure. 



We want to represent the fermion integral of eq. ([4.8| ) as an integral with the effective 
fermion action, eq. (|3.35|) . The effective action was derived from a product of /'s and 
N_ J's. In eq. (|4.8|), however, only N± — 1 of the /(/)'s are "free", while one /(/) is 
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connected to the operator, J-'[A]. The effective operator in the sense of eg. ([4.5|) can thus 
not be immediately identified with the vertices VV±) eq-( [4.11|) . Rather, it is of the form 



(4.12) 



where TZ±[ilj^,ip] are functions of the fermion field, which compensate the fact that we 
have "taken out" an /(/) from the N± /(/)'s that make the effective action, eq. (|3.15 ). 



To determine the factors Tl±, we proceed with the evaluation of eq. (|4.8|) in analogy to 
the derivation of the effective action, eq.( |3.15D . Using the same integral representation for 
the powers of one-instanton vertices, eq.( p.21| ), but this time for N± — 1 "free" vertices, it 
is easy to see that TZ± must be chosen as times a factor which compensates the shift in 
the saddle point compared to the integral defining the effective fermion action, eq.( |3.2lD . 
After a straightforward calculation one obtains 
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iV2 



F:p exp 



Y 



(v;>ei,{v-) 




1 - 



(Y-) 



eff. 



(4.13) 
(4.14) 



Here, Kj-f^/j^, ^/'j are the 't Hooft vertices of the effective fermion action, eq. ( |3.36|) . In the 
second equation we have used that, to leading order in l/N^, the vacuum average of Y± is 



eff 



(4.15) 



as follows from a straightforward calculation, using the saddle point equation, eq. (|3.25| ). 

The effective operator, as defined by eqs.( [4.11"| , [4.12| , ^.14| ), has an intuitively obvious 
form. It consists of 't Hooft-type vertices describing the interaction of the gluon operator 
with the fermions mediated by one instanton. The interaction strength is governed by 
M, the dynamical fermion mass. The one-instanton vertices are multiplied by certain 
universal functions of the fermion fields, TZ±. The meaning of these functions can be 
elucidated by expanding eq. ( [4.14| ) around the saddle point. To first order in Y± — {Y±)es, 
one has 



7^^ 



1 



1 



y± - {y±) 



eff 



eff 



{y±) 



eff 



y-1 



(4.16) 



i.e., the functions 1Z± can be regarded as the "inverse" of a 't Hooft vertex in the vicinity 
of the saddle point. As will be seen below, the presence of the functions TZ±[%p\il)\ in 
eq. ([4.12|) , in addition to normalizing the effective operator, has important consequences 
when considering the connected average of the effective operator with other fermionic 
operators, such as hadronic currents. 

Let us calculate the average of the effective operator, eq. ( [4. 21] ), in the vacuum of the 
effective fermion theory. To leading order in l/N^. only the disconnected average of the 
factors of eq. ([4.12|) needs to be considered. Furthermore, when averaging the exponential 
factor of eq. ( ^4.14|) we may use 



(X") 



eff 



(^)e"f 



eff 5 



(exp X) eff = exp (X) eff 



(4.17) 
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{X is an arbitrary function of the fermion fields). With eq.( |4.15| ), we obtain 



and thus 



eff 



eff 



eff 



(Y-) 



eff 



(4.18) 



(4.19) 



As expected, the average of the effective operator is given by the normahzed average of 
the one-instanton vertices, multiphed by the number of /'s and /'s in the ensemble. 

Let us construct the effective operators corresponding to the simplest gauge-invariant 
gluonic operators, F^^ and F^^F^y. For one /(/), 



192p4 



n4 ' 



[p^ + {x — 

and the corresponding one-instanton vertices, eg. ( H^.llQ , are 



F^uF^y{x) — ihi^ (x) 



(4.20) 



Yf^±{x) 



Ypp jJyX) 



IV 



Nf~l 



N J 



192p^ 



[p2 + (x - z) 



det J±{z) 



(4.21) 
(4.22) 



In practice, we are interested only in matrix elements of these operators at zero momentum 
transfer (forward scattering), and it is sufficient to consider the integral of the local 
operators over the four-volume. The integrals of eqs.( [4.21"] , [4.22| ) are identical to the 
't Hooft vertices, eg. (|3.36|), times the action of one instanton. 



327r2y±, 



^ Ypp jJyX) 



±327r2r±. 



The vacuum averages of the effective operators, eq.( [4.19| ), are thus 



eff 



327r2(A^+ + 



(TxF, 



flU 
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fixed -Ar± 



(4.23) 



(4.24) 



and similarly for F^j^F^^,^ with + replaced by — A^_. The average of the 
effective operator in the vacuum of the effective fermion theory is thus identical to the 
average of the original operator in the fixed-A^± instanton ensemble, in accordance with 
the definition, eg. ([4.5|). 

Consider now vacuum averages of the effective operators with hadronic currents. An 
important property of the effective operator / d^x "F^^" is that its average together with 
an arbitrary fermionic operator, Ol^p^ip], reduces to a disconnected average. 



o I d^x 



eff 



((^)eff 

327r2(A^+ + N^){0) 



eff 



(4.25) 
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The connected part of the average, which would be of the same order in l/N^ as the 1/Nc- 
correction to the disconnected part in eg . (14.251 ), is zero. (Here, the terms "connected" and 



'disconnected" refer to the effective fermion theory.) To see this, let us calculate explicitly 



the connected part of eq.( [4.25|) . The operator O can connect to any of the factors in the 
effective operator, eg. ( [4.211 ). The contribution of O connecting to the functions TZ±, 



eg. ([4.14|) , is calculated using the rules 

(0,X")eff,conn = n{0 , X) ,s,conn {X):^\ (4-26) 
(C',expX)efr,conn = (C, X)efr, conn f ^ exp X j = (O, X)efr, conn exp (X)eff, 

which apply to the leading order in l/N^. One easily obtains 

((^,>^±^±)eff,conn = (^^^ >^±)efr, conn (^±)efr + (>^±)eff (^^^ ^±)efr, conn = 0,(4.27) 

in accordance with the interpretation of TZ± as the "inverse" of Y±, eg.( [4.16| ). From this 
it follows that 



{O J rf^x"F2;')eff,conn = 0. (4.28) 



An analogous property holds for connected averages with the operator / d'^x '^F^^F^^" . 

The role of the effective operators / d'^x "-F^^" and / d'^x ''F^^F^y" , when averaged 
together with other fermionic operators, is simply to multiply the average of these oper- 
ators by 327r^(A^+ ± A^_). There are no other dynamical effects from the insertion of an 
operator / d'^x "-F^j," or / d'^x ^^F^yF^^^ in a correlation function in the effective fermion 
theory. This is what one should expect: the original operators / d'^xF^^ and / d'^x Ffj_yF^y 
measure just the number of J's and J's in the ensemble, and, by definition, the average 
of the effective operators in the effective theory is egual to the average of the original 
operator in the instanton ensemble with fermions. 

The absence of a connected average in eg. (|4.25|) has important conseguences when 



passing from canonical to grand canonical averages. This property is crucial for realiz- 
ing the trace and f/(l)yi-anomaly in the matrix elements of the operators / d^x F^^ and 
/ d^xFfj^yF^y in hadronic states, as will be seen in section 5. 

4.2 Fluctuations of the numbers of instantons 

The effective fermion action, together with the effective operators for gluons introduced in 
the previous section, allows to calculate normalized averages of operators in the ensemble 
with a fixed number of instantons. Physical guantities, however, are given by averages over 
the grand canonical averages, described by the grand partition function with fermions, 
eg.(|2.1|). To pass from canonical to grand canonical averages, we have to sum the fixed- 



N± normalized average, egs.( [4.2| , [4.4|) , over the distribution of the numbers of J's and J's 
in the ensemble. The weight for a configuration with given A^^, A^„ is given by the value 
of the fixed-A^± partition function with fermions. With the variational approximation, 
eg. (13), 

P(Ar+,Ar_) oc = Z^v^Dit^^. (4.29) 
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(For convenience, we shall assume the distribution P{N+, N^) to be normalized.) Since 
we consider only quadratic (gaussian) fluctuations around the equilibrium values, = 
= ^{N), we may instead of N± equivalently average over independent distributions 
oiN = N+ + N_ and A = N+ - N_. 

The normalized grand canonical average of a fermionic operator, Olip^ip], is defined 

as 

{0[^/j\^]) ^ Yl P{N^,N^) {0[ij\ij]) (4.30) 

N+,N. ° 

For a gluon operator, J^[v4], the only difference is that the canonical average is calculated 
with the help of the corresponding effective operator, 

{T[A]) ^ 5: P(iV+,iV_) ("^"[^t^^]\ (4.31) 

Here, it is assumed that the fixed-A^± averages over the effective fermion theory have been 
evaluated as functions of N±, in the vicinity of = = \{N) and to first order in 
A = N+-N_. 

In particular, with this definition, the grand canonical average of / d'^xF^^, using 



cfT 



eq.(pl), is 

= J2 P(iV+,iV_)(iV+ + iV_) = {N). (4.32) 

N+,N. 

The effective operator technique reproduces the usual value of the vacuum gluon conden- 
sate. In the same way one obtains that the grand canonical average of / d'^x F^j^yF^^ is 
zero. 

The relative width of fluctuations of — (A^) and A is inversely proportional to the 
size of the system, cf. eqs.( |2.20| , |3.56| ). Thus, fluctuations have no effect on averages of 
intensive quantities in the thermodynamic limit. For example, the correlation function 
of two baryon currents in the grand canonical ensemble simply reduces to the one in the 
canonical ensemble, eq.( [4.2|) , evaluated at the equilibrium value of iV±, 

{Mx.VUx,)) = (^^(-i)4(a:.))^,|,^^,,,^„ + (4.33) 

As a consequence, the infrared-singular term of the effective fermion action proportional 
to A, eq.( |3.35| ), does not affect free hadron correlation correlation functions, as it should 
be. However, fluctuations of A^ and A lead to non-trivial consequences, when considering 
connected grand canonical averages of fermionic currents with extensive operators, such 
as / d'^x F^y and / d^x F^yF^y. We shall now discuss these effects, considering separately 
non-topological (A^ — (A^)) and topological (A) fluctuations. 

Non-topological fluctuations. Let us consider the connected grand canonical average 
of an arbitrary fermionic operator, 0[ip\tp], with the operator / d^x F^^. According to 
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2 J5 



Cffy 



(It is understood that A = 0; we do not write the averaging over A here.) Here, we have 



employed the property eq.(|4.25|) for the fixed- average. In the last equation, we have 



used the fact that the relative width of the fiuctuations, {{N — (N))"^) / {N), becomes small 
in the thermodynamic limit. (In addition, the width is proportional to 1/Nc, eq.( 2.2U| )). 



In eq. ( |4.34|) it is implied that {0)es is first evaluated by integrating over the fermion 
fields for fixed N, and then differentiated with respect to N. This differentiation also 
includes the dependence of {0)es on through the instanton size, p, which is related to 
N/V according to eq.( p.44|) . Since the only scale in the effective fermion theory is the 
instanton density, N/V, the operator N{d/dN) measures the total scale dependence of 
{0)es- In this sense, it plays the role of A{d/dA) in QCD. 

Taking the A^-distribution from gluodynamics (quenched approximation), eq.( |2.20| ), 
eq. ([4.34|) becomes 

We have reduced the connected grand canonical correlation function to the A^-derivative 
of a fixed-A^ correlation function, calculable in the effective fermion theory at equilib- 
rium, A^ = (A^). Eq.( 4.35| ) is the statement of the fact that the renormalization (scale 



dependence) properties of QCD are reahzed in our approach. As a special case, eq.( |4.35| ) 
contains the low-energy theorem of scale invariance for the vacuum averages of local 
fermionic operators If (9 is a local operator, then, as N/V is the only scale in the 

effective fermion theory. 

Here, d is the naive dimension of the operator. 

We emphasize that eq. (|4.35| ) is valid also for non-local fermionic operators, in partic- 



ular, for operators of the form O = Jn{.Xi)J'n{.X2)- Such operator have an intrinsic scale, 
defined by the distance of the two points, xi,X2, and the vacuum average is no longer 
a homogeneous function of N/V. However, eq.( [4.35D still leads to useful results in the 



limit of large distance, \xi — X2I — > 00. In this regime, the correlation function possesses 
an asymptotic expansion, the coefficients of which are dimensionful and therefore propor- 
tional to powers of the instanton density. This fact can be used to derive the nucleon 
matrix element of / d^x F"^^ (see section 5). 
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To summarize, a correct realization of the scale dependence properties of QCD is 
achieved in our approach already in quenched approximation, neglecting the influence 
of the fermion determinant on the distribution of A^. This is different for topological 
fluctuations, which we shall now discuss. 

Topological fluctuations. We now consider correlation functions with the topological 
charge operator, / d'^x F^^Ffj^^,. In this case, the average over fluctuations of A = A^-|. — 
leads to non-trivial consequences. In analogy to eq.(|4.35|), the grand canonical average of 



an arbitrary fermionic operator, 0[ip\tp], with the topological charge becomes 

= i;nA)A{OU = (A=> (OU) . (4.37) 

Here, we have used that (A^)^/^/A^ is small for large A^. The derivative with respect to 
A of the fixed- average in the effective theory can be calculated as 

(4^(0U] = {O^-^} (4.38) 



dA 



A=0 



dScs I 




[^^^] = W77 E-7M(^+-^-)- (4-39) 



dA ' N,VM 



Here, dSes / dA[ip\ ip] is treated as an operator insertion, and the average over the effective 
fermion is now to be taken with A = 0. 

For the dispersion of A in eq.( |4.37D , we must now use the distribution which includes 



the effect of the fermion determinant, eq. ( p.55[ ). With this dispersion, one obtains 



O^^jd'x F,,F,}j = {O (n - r_)),g I A=o. (4.40) 

Notice that the 1 /m-divergence in the chiral limit, which is present in the effective fermion 
action for A 7^ 0, is compensated by the fact that the dispersion of A, eq.( p.55|) , is of 



order m. As a result, the right-hand side of eq. ( 4.401) is finite in the chiral limit. It would 



be senseless to average over A using the quenched distribution, (A^) = (A^) — the result 
would be divergent in the chiral limit. This is an important difference to the description 
of non-topological fluctuations. 

We now want to demonstrate that eq. ([4.40|) is consistent with the U{1)a anomaly of 
QCD. To this end, let us derive the f/(l)-axial current and its divergence in the effective 
fermion theory defined by eq.( p.35|) , at A = 0. Consider a [/(l)-chiral rotation with 
parameter e, 

— i> exp (z£:75) '?/', if)'^ — >• ■^/'^ exp (^£75) . (4-41) 

The change of the effective action, eq. (|3.35|) , gives the divergence of the axial current, 

jd'xd,.h,{x) = ^1^^^ = 2Nf{Y^-Y_). (4.42) 
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The current itself is obtained from the fermion kinetic term, as usual^ 



= E^/(^)7/.75V'/(a;). (4.43) 
/ 



Thus, eq.(l^) can equivalently be written as 



O 



3271- 



jd'xF,^F,}) = lo-l- j d'xdMx)\ 1^^^. (4.44) 

\ f I off 



This is precisely the relation which in QCD would follow from the anomaly equation, 
eq. (|3.38|) , in the chiral limit. We have thus shown that our prescription, based on the 



effective fermion theory and the effective operator for / d'^x F^^F^^, fully realizes the U{1)a 
anomaly at the level of correlation functions. 

One should note the different character of the U{1)a symmetry breaking in QCD and 
the effective fermion theory. In QCD, the symmetry is broken anomalously, i.e., by the 
high-momentum components of the fermion field. The effective fermion theory describes 
only the low-momentum components of the fermion field, up to a cutoff, and so, by 
definition, it has no anomaly. (We have assumed this in deriving the divergence of the 
axial current, eq.( [4.42D .) However, in the effective theory the U{1)a symmetry is broken 
explicitly, which results in a non-zero divergence of the axial current. 



5 Hadronic matrix elements of and F^jjF^jj 

We now apply the techniques developed in the previous sections to the study of the nucleon 
matrix elements of the operators and F^^F^y. On general grounds, the nucleon matrix 
elements of these operators can be parametrized as 

{p'\FlM\p) = As{q^)mr,u'u, (5.1) 
{p'\F^uF^,{{))\p) = Ap(g2)m^u'z75U, (5.2) 

where q = p' —p, and u, u' denote nucleon spinors. In QCD, ^5(0) is uniquely determined 
by the conformal anomaly ||34|| , 



M^) = ^, (5.3) 

whereas Ap(0) is related to the isosinglet axial coupling constant, g^^ ^ by virtue of the 
axial anomaly equation, eq.( p.38|) , 

Ap(0) = (5.4) 

{p\Jb^p) = gA^ui^lbU. (5.5) 



^We neglect the contributions to the current due to the momentum dependence of the dynamical 
fermion mass. They are parametrically small, of order Mp. 
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We show now that these relations are reahzed in our approximate description. Thus, the 
renormahzation properties of QCD and the axial anomaly are correctly taken into account 
even at the level of hadronic matrix elements. 

Consider first the operator F^^{x). The nucleon matrix element in our approach is 
extracted from the large-time limit of the connected euclidean correlation function, 

MT)jU-T)FUx)) - (jiv(T)4(-T)) (f'(x)) . (5.6) 



(Here, "connected" refers to the full grand canonical average.) For brevity, Jn{T) = 
Ji\f{0,T), etc. Since we are interested in the zero momentum transfer matrix element, we 
integrate the operator F^^{x) over the euclidean four-volume. One obtains 



As{0) = lim 



MT)4{-T) J d'xFl) - {Jn{T)JI{-T))Q 



d'xF% 



T^oo 2m^T (j^(T)4(-T); 

(5.7) 

By virtue of the effective version of the "low-energy theorem", eg. ( [4. 351) , this becomes 

^5(0) = hm l^327r^| fivAiog(j^(T)4(-T)) ) . (5.8) 
T-oo 2mNT h \ dN \ ' ) n={n) 

For large T, the euclidean correlation function decays exponentially, 

{Jn{T)4{-T))^^ ~ exp(-2Tm;v). (5.9) 

Here, tjin is the nucleon mass, as calculated in the effective fermion theory for a fixed 
number of instantons, N . For dimensional reasons, 

oc ( — I , (5.10) 

since the instanton density is the only scale in the effective fermion theory. (The instanton 
size, p, is related to N jV by eq. (|2.44|) .) With eg. ( |5.9| ) and eg. (|5.10|) , we obtain from 
eg. (1:1) 

397r2 

A.(0) = (5.11) 

in agreement with the trace anomaly. This result is again a manifestation of the fact that 
the instanton vacuum preserves the renormahzation properties of QCD. 

Let us now consider the nucleon matrix element of the topological charge, eg. (|5.2|) . In 
analogy to eg.( [5.7| ), it is obtained as 



. ^Jn{T)JI{-T)F,^F,Ax), 
Ap(0) = hm — ^ -. T ^ '-. (5.12) 

^ ^ 2m^^ (^iv(r)4(-r)) ^ ^ 
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With the reahzation of the axial anomaly in the effective theory, eq.( [4.44| ), we have 

167r2 

1 

Ap{0) = lim 



Nf 



cflF 



A=0 



MT)Jli-T) 



cS 



(5.13) 



A=0 



The correlation functions here are to be evaluated in the large- Ai'c limit, in which the 
nucleon is described as a chiral soliton in the effective fermion theory |2ll]. In momentum 
representation, the free nucleon correlation function develops a pole at m^v, and a nucleon 
state emerges in the effective theory. Saturating the correlation functions numerator and 
denominator of eq. (|5.13| ) with one-nucleon states of the effective theory, one finally obtains 



Ap{0) 



327r2 (0) 



(5.14) 



where g^^^ now denotes the axial coupling of the nucleon, as evaluated in the effective 
fermion theory from the matrix element of the isosinglet axial current, cf. eq.( ^.5|) . We 
have thus shown, that the axial anomaly is realized in our effective description even at 
the level of hadronic matrix elements: the nucleon matrix element of F^yF^y reduces to 
the nucleon isosinglet axial coupling within our effective description. 

The nucleon isosinglet axial coupling, g^^ , has been evaluated numerically in the chiral 
quark soliton model of the nucleon, which is obtained from the above effective fermion 
theory by bosonization Calculations give g*^^ ~ 0.36 which agrees well with the 



value obtained in a recent analysis of polarized nucleon structure functions^] [p5|] . 

It has been suggested by Anselm ||3^ that a small value of g'^^ could be explained 



by equating Ap{Q) with ^5(0), and thus g'^^ with —Nf/b, at the constituent quark level. 
He argued that Ap{0) = ^5(0) could be obtained from the fact that F^^F^^ = ±F^^, 
for one /(/). Our derivation shows that there is no justification for such a picture. 
The fluctuations of the numbers of instantons in the ensemble are crucial to consistently 
describe the matrix elements of F^^ and F^yF^y. For a single constituent quark instead 

of the nucleon, we immediately obtain from eq.( p.l3| ) g^^ = 1. The fact that a value of 



g^^^ considerably smaller than 1 is obtained for the nucleon 



can be understood as an 



effect of the meson cloud of the nucleon (polarization of the Dirac sea by the Nc valence 
quarks) . 



6 Conclusions and outlook 

In this paper we have considered a description of non-perturbative phenomena in QCD in 
the form of a grand canonical ensemble of instantons. The fluctuations of the number of 
pseudoparticles play a crucial role in realizing the renormalization properties of QCD in 
the resulting statistical mechanics system. Fluctuations of = A^++A^_ (non-topological) 
are related to the trace anomaly, while fluctuations of A = — (topological) realize 
the t/(l)^-anomaly. The relation between the QCD anomalies and the fluctuations of N± 

^For a review of the "proton spin problem" , see p4[ . 
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has two directions: First, starting from the QCD low-energy theorem of scale invariance 
and the anomalous U{1)a Ward identities, one can predict the dispersion of N and A 
in the instanton medium from first principles Second, as we have seen, the instanton 
medium, obtained with a concrete ansatz for the instanton-instanton and instanton- 
fermion interaction, exhibits a dispersion of and A in accordance with the general 
principles. An essential requirement is a fully self-consistent treatment of the instanton 
medium, based on instanton interactions — any external stabilizing mechanism would 
violate the renormalization properties. 

Our second aim has been to develop a consistent practical prescription to evaluate op- 
erators — in particular, gluon operators — in the grand canonical ensemble of instantons 
with fermions. We have invoked the large- A^^c limit, which simplifies the dynamics in two 
important aspects: First, the relative width of the fluctuations of the total number of 
pseudoparticles becomes small. Second, in the large-A^c limit one can explicitly construct 
the effective fermion action for an ensemble with a fixed number of instantons. This en- 
ables us to evaluate correlation functions in a two-step procedure: We first integrate over 
instanton coordinates and the fermion field in a canonical ensemble, and then average 
the result over fluctuations of N±. All fixed- correlation functions can be represented 
as integrals over the effective fermion theory, in which gluon operators are replaced by 
effective fermion operators. We have shown that in this approach the trace and U{1)a 
anomalies are realized at the level of hadronic matrix elements of the operators / F^^^ and 
/ F^tyF^i,. Consequently, the result for the isosinglet axial coupling of the nucleon of [p^ 
may be regarded as a consistent estimate of the nucleon matrix element of the topological 
charge in the instanton vacuum. 

We have established, that our scheme of approximations (variational treatment of 
instanton interactions, A^^c-limit for the fermions, zero mode approximation) provides a 
framework for evaluating gluonic operators, which is consistent with the basic renormal- 
ization properties of QCD. Any attempt to improve these approximations {e.g. using 
more sophisticated trial functions, including explicit correlations between instantons due 
to fermions, or going beyond the zero-mode approximation for the fermion propagator) 
should preserve these standards of consistency. 

An application of the methods developed here is the calculation of nucleon structure 
functions, both for leading and non- leading twist. By the operator product expansion of 
QCD, the moments of the structure functions are related to forward matrix elements of 
certain local operators which contain gluon fields. These matrix elements can be evalu- 
ated using the technique of effective operators developed here. The moments are obtained 
at the scale set by the instanton background, typically 1/p, and have to be put into QCD 
evolution equations to compare with experimental information at higher momenta. Work 
in this direction is in progress. 
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